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Abstract 

A family of nonlinear ordinary differential equations with arbitrary order is 
obtained by using nonextensive concepts related to the Tsallis entropy. Appli- 
cations of these equations are given here. In particular, a connection between 
Tsallis entropy and the one-dimensional correlated anomalous diffusion equa- 
tion is established. It is also developed explicitly a WKB-like method for 
second order equations and it is applied to solve approximately a class of 
equations that contains as a special case the Thomas-Fermi equation for an 
atom. It is expected that the present ideas can be useful in the discussion of 
other nonlinear contexts. 
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I. INTRODUCTION 



As it is well known, the equation 

dp 



whose general solution is 



p{x) = Poexp(-Ax) , (2) 



has many applications. In general, Eq. (|T|) models systems related to an extensive context 
(in the thermodynamical sense). Typical examples are the systems based on independent 
events, for instance, the radioactive decay of noninteracting nucleus. In this context, it is 
interesting to observe that the solution of Eq. (|l|) has an entropic interpretation. In fact, if 
the usual entropy, 

rb 

S = — p{x) Inp(x) dx , (3) 

Ja 

is maximized subject to the constraints 

rb 

a = p{x) dx (4) 



and 

rb 

P = X p{x) dx , (5) 

J a 

the solution of Eq. (|l]) for a < x < 6 is obtained again. The parameters a and (3 are 
adjusted in order to obtain po and A. If a = 1, p{x) can be interpreted as a probability. 
However, this choice for a is not necessary because all the conclusions presented here are 
independent of the a value. 

Now, it is possible to ask how Eq. ([1|) must be generalized in order to describe depen- 
dent events {nonextensive context). As a guide to answer this question it will be considered 
in this work a generalization of the entropy (|^) employed by Tsallis (see also Ref. 0) 
in a nonextensive statistical mechanics. The above choice is motivated by the fact that 



this generalized entropy (Tsallis entropy) has been apphed successfuly in the discussion of 
many situations where the nonextensivity plays an important role, for instance, Levy su- 
perdiffusion and correlated anomalous diffusion , turbulence in two-dimensional pure 
electron plasma [0], dynamic linear response theory ^ and Green functions 0, perturbation 
and variation methods for calculation of thermodynamic quantities |TD[ , low-dimensional dis- 



sipative systems |]TT|, simulated annealing and optimization techniques [|12| and connection 



with quantum uncertainty [|l3l. In other words, the main purpose of this work is to obtain 
generalizations of Eq. (|l|) based on nonextensive Tsallis entropy, giving physical applica- 
tions. More specifically, it is obtained a generalization of Eq. ([l|) by using Tsallis entropy 
(Section II), and applications to the motion of a particle in a fluid medium and chemi- 
cal kinetics are presented. A further generalization for higher order nonlinear differential 
equations is given (Section III) and it is applied to discuss a new connection between the 
correlated anomalous diffusion equation and the Tsallis entropy. Moreover, it is introduced 
a WKB-like approximation in order to study a family of nonlinear differential equations 
based on Tsallis entropy (Section IV), and finally this approximated method is employed 
to the case of the Thomas-Fermi equation. 

II. FIRST ORDER NONLINEAR DIFFERENTIAL EQUATION 

The Tsallis entropy can be written as 

J a 1 — (J 

Furthermore, the constraint (^ is currently substituted 0,0] by 

(3= [\p{xydx (7) 



and the constraint (1) remains unchanged. When the entropy (|]) is maximized subject to 
the constraints @) and (|^, p{x) can be written as 

p{x) = Poll - {I - q) pI'^ Xx] . (8) 



This distribution generalizes the exponential one (see Eq. (|D). In the above expressions the 
parameter q E TZ characterizes the degree of nonextensivity, in particular the entropy 
and the constraint (^) are obtained as limiting case when q —>■ 1, recovering the extensive 
case. A remarkable fact is that Eq. (^) can be applied directly to several contexts successfuly. 
In addition to other works already cited, the study of solar neutrinos and Zipf law 
are relevant examples of such applications. 

By direct inspection it is verified that the above generalization of the exponential satisfies 
the nonlinear equation 

| + A,' = (9) 

subject to the initial condition p(0) = Po- This equation is the generalization of Eq. (|Tp 
based on the nonextensive Tsallis entropy. When the notation Ae// = Xp'^~^ is employed, 
Eq. can be interpreted as a decay with memory, thus q is in some sense related with the 
memory of the system. 

A. Applications 

Before generalizing Eq. (Q) to higher orders it is illustrative to present applications for 
it. One of them describes the mean motion of a particle in a fluid medium without external 
force. In this case, the motion equation can be written as 

mdv/dt = -hv\ (10) 

where v and h are respectively the velocity and the friction coefficient of the particle relative 
to the medium and m is the particle mass. Here, the parameter q is related to the turbulent 
flow (nonextensive behavior) and for q = 1 (slow motion) the extensive behavior is recovered. 
Another example comes from chemical kinetics. In this case, the concentration Ca of a given 



species A obeys the empirical equation |T6| 



dCA/dt = KC'XC^^Cl--, (11) 



where K is the reaction constant and a, /3, 7, • ■ ■, refer to the concentration of chemical 
species A, 5, C, ■ ■ ■, present in the reaction. In this expression, a, /3, 7, ■ ■ ■, are respectively 
the order of the reaction with respect to A, 5, C, ■ ■ ■, and the sum a + /5 + 7 + ■ ■ ■ is the 
overall order of the reaction. In some cases the concentrations Cb-, Cc, ■■■ can be considered 
constant, thus the above equation reduces to the form of Eq. (|^) with q = a. In this way, 
the parameter q becomes the order of the reaction for the species A. 



III. FAMILY OF NONLINEAR DIFFERENTIAL EQUATIONS 

The natural generalization of Eq. (|1|) for higher order differential equations, in the sense 
that p = Poexp(— Ax) is a particular solution, is the family of linear differential equations of 
arbitrary order with constant coefficients. The iV-order element of this family is 

^ d"n 

E«n^ = (12) 

and the corresponding algebraic equation for A is 

N 

^ a„(-A)" = . (13) 

ra=0 

In view of the previous remarks, the generalization of these ideas in the Tsallis context is 
based on the replacement of the particular solution p{x) = poexp(— Ax) by p(x) = Po[l — 
(l-g)p^-Ux]i/(i-^). Thus, 

n=0 

represents the nonlinear ordinary differential equation of constant coefficients that general- 
izes Eq. ([1^) and the corresponding generalization of Eq. (|^) is 



N-2 JV-1 N-2 

aNi-Xf n iiJ + 1)? - + E «n(-A)" n KJ + 1)? - j] + ao = . (15) 

j=0 n=l j=N-n-l 

This equation is applicable for N > 2 and when = 1 Eq.(]T3|) must be replaced by 
ai(— A) + ao = 0. Without loss of generality the A^-th term of Eq. (|Tip, d^p/dx^, was 



considered linear in p{x). As expected, Eqs. (|Tj) and (^) reduce respectively to Eqs. (|T2D 



and (|I3D in the limit g — 1. Furthermore, it is important to remark that a superposition of 
particular solutions of Eq. (|14|) is not another solution of this equation, in contrast with the 
case of Eq. (0). For instance, the generalization of circular and hyperbolic functions based 
on the superposition of Eq. (||) ||T^ is not a solution of Eq. Thus, in the following 

applications only particular solutions previously described will be considered. 



A. Application to correlated anomalous diffusion 

It was recently shown that the one- dimensional correlated anomalous diffusion equation 
(generalized nonlinear Fokker-Planck equation without external force), 

has a solution related to the Tsallis entropy 

This solution is a generalization of the Gaussian one, reducing to it in the limit q — >■ 
1. On the other hand, by using the generalization developed above for linear ordinary 
differential equations, a new relation between Tsallis entropy and correlated anomalous 
diffusion equation can be obtained. 

The new connection is based on the substitution of 

(Pix,t)=Tit)X{x) (17) 

into Eq. (|T6D. This separation of variables leads to dT/dt = —a T'^ and d^Y/dx'^ = 
—a Y^^" , where a is a separation variable constant and Y = X'^ . Since these equations are 
members of the family ruled by Eq. (0), a new connection between the one- dimensional 
correlated anomalous diffusion equation and the Tsallis entropy is established. Furthermore, 
as a consequence of this general procedure, a set of particular solutions of the anomalous 
correlated diffusion equations is obtained, namely. 



x,t)=ro l-{l-v)T^-^ at Xo 



1 - u 



1 + ( ^ ) Xt-'^h X 



2/(^-1) 



(18) 
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where Tq and Xq are constants and A is the solution of the equation (1 + z/)A^ + 2z/cr = 0. 
Note that the above procedure leads to complex solutions for a > and it is a natural 
extension of the method of separation of variables applied to the usual diffusion equation 
{u = l). 

The previous method can be easily applied to the generalized diffusion equation recently 
proposed by Tsallis and Bukman f^, 

-W = ^ (19) 

In fact, it is sufficient to replace (p^ by ip and u by z//yU. In general, the connection between 
the porous media diffusion equation and Tsallis entropy, based on separation of variables 
and the family of nonlinear partial differential equation (Eq. (p!4D), can be easily extended 
for other nonlinear partial differential equations. 

IV. NONLINEAR DIFFERENTIAL EQUATIONS AND WKB-LIKE 

APPROXIMATION 

Eq. @ can be generalized further if we allow A to become a function of x. In this case, 
the solution of Eq. (H) with p{0) = po becomes 

1/(1-9) 



p{x) = Po 



l-(l-g)prW A(^)d^ . (20) 

Jo 



In a similar way, we can allow that the constants a„ in Eq. ([I^) become functions of x. This 
procedure leads to 

N jn 

Ean(x);f^p(^-")(-^)+^ = 0. (21) 
n=0 '^^ 

As in the linear case (g = 1) there are no general solutions for these equations. Consequently, 
it is natural to perform approximated analyses to obtain some information about the solu- 
tions of Eq. (pTD . In the following discussions, among other possibilities, a generalization of 
WKB method for g 7^ 1 is developed explicitly for the N = 2 case. 
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As it is well known, in the WKB method [|T8[ approximate solutions of equation 
d^p/da;^ = f{x)p can be obtained when f{x) is a slowly varying function. In this case, 
it is employed an auxiliary function g{x) defined by the relation p{x) = exp {g{x)). /^From 
the previous developments, a natural generalization of these equations are respectively 



dx"^ 



f{x)p 



(22) 



and 



p(x) = + . 
In terms of g{x), Eq. ( ]22|) becomes 



(23) 



(24) 



Following again the usual WKB approach, the term with d^g/dx'^ is neglected in the first 
approximation, hence 



g{x) = ±q-^'^ j f{xf'^dx. 
Consistently, the validity condition of this approximation is 



(25) 







d//dx 


< 


/ 


dx'^ 




2gV2/i/2 


l + il-q)g 



(26) 

Since Eq. (|22|) is nonlinear for g 7^ 1, the superposition of solutions is not a solution, therefore 
the present development is indicated for situations where some particular solutions can be 
considered as good approximations. Notice also that the solutions (p5|) can be improved 
through iterations (replacing successively improved Eq. ( ^5]) into Eq. (p4D). 



A. Application to Thomas-Fermi equation 

To exemplify the previous development, we consider the Thomas- Fermi equation for a 
free atom [|r^ 
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2/3/2 _ (27) 



In the free neutral atom case the boundary conditions are y{0) = 1 and y{oo) = 0. By 
comparing Eq. (13) with Eq. (|2|) it is verified that /(x) = x"^/^ ^nd g = 5/4. Choosing 
the particular solution adjustable to these conditions it becomes 

^3/4^ 

In this example, the validity condition (EB) becomes 



y{x) =(l + ^ ) . (28) 



1 + M a;-3/4 15 _ (29) 



This condition indicates that the approximation is better for larger x. Furthermore, Eq. 
is in satisfactory agreement with a numerical calculation (see, for instance, Ref. ||19|). In a 
general context, when f{x) is a smooth function, the corresponding approximate solution 
becomes more accurate. In particular, when f{x) is a constant this solution becomes exact. 



V. CONCLUSIONS 

Summing up, the nonextensive concepts based on the Tsallis entropy were employed to 
obtain a family of nonlinear ordinary differential equations. The first order equation of this 
family is a nonextensive generalization of the exponential decay equation. Moreover, by 
using a separation of variables procedure and the above family of equations, a connection 
between the correlated anomalous diffusion equation and the Tsallis entropy is obtained. In 
addition to this, for second order equations we presented a WKB-like approach to obtain 
approximated solutions. This procedure was used in the context of the Thomas-Fermi 
equation for a free neutral atom, and it was shown that the well known solution precisely 
correspond to g = 5/4. In general, the developments introduced in this work indicate 
that many nonlinear effects are closely related with nonextensive concepts in the Tsallis 
framework. 
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